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Abstract 

Thermodynamics provides powerful constraints on physical and chemical systems in 
equilibrium. However, non-equilibrium dynamics depends explicitly on microscopic 
properties, requiring an understanding beyond thermodynamics. Remarkably, in di¬ 
lute gases, a set of universal relations is known to connect thermodynamics directly 
with microscopic properties. So far, these “contact” relations have been established 
only for interactions with s-wave symmetry, i.e., without relative angular momentum. 
We report measurements of two new physical quantities, the “p-wave contacts”, and 
present evidence that they encode the universal aspects of p -wave interactions through 
recently proposed relations. Our experiments use an ultracold Fermi gas of 40 K, in 
which s-wave interactions are suppressed by polarising the sample, while p-wave in¬ 
teractions are enhanced by working near a scattering resonance. Using time-resolved 
spectroscopy, we study how correlations in the system develop after “quenching” the 
atoms into an interacting state. Combining quasi-steady-state measurements with 
new contact relations, we infer an attractive p-wave interaction energy as large as half 
the Fermi energy. Our results reveal new ways to understand and characterise the 
properties of a resonant p-wave quantum gas. 
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A fundamental question provoked by observation of natural systems is how macroscopic 
and collective properties depend on microscopic few-body interactions. Ultracold neutral 
atoms provide a model system in which to explore this question, since in certain conditions, 
few-body interactions can be tuned and characterised precisely. Over the last decade, a direct 
link has been made in these systems between thermodynamic properties and the underlying 
isotropic (s-wave) interactions. At the centre stage is a quantity called the “contact” P- 
[6], which describes how the energy of a system changes when the interaction strength is 
changed. Surprisingly, the contact is also the pivot of a set of universal relations, that 
constrain numerous microscopic properties, including the two-particle correlation function 
at short range. These relations apply regardless of temperature, density, or interaction 
strength IM E], to fermions and bosons and in one-, two-, and three-dimensional 

systems [7[ 1101412] . Contact relations have also been extended to Coulomb gases [13] and 
neutron-proton interactions [lTj. Despite the breadth of this discussion, measurements of 
the contact have so far been restricted to systems with s-wave interactions. 

In general, the relative wave function of any pair of particles can be decomposed into 
components with angular momentum equal to an integer multiple £ of h quanta. In a spin- 
polarised Fermi gas, quantum statistics forbids short-range interactions with even values of 
l. Therefore, the first allowed scattering channel has £ = 1 (p-wave), which is typically weak 
due to the centrifugal barrier (see Fig. [I]): the scattering cross section decreases with the 
square of the collision energy eg. However, resonant enhancement of p -wave collisions has 
been observed in 40 K and 6 Li [T61IT8] . raising the possibility of studying a gas with strong, 
tuneable p-wave interactions. 

Interest in p -wave systems originated with liquid 3 He, which at low temperature is a 
superfluid of pairs with £ = 1 in a spin triplet state. It is thought that Sr 2 Ru0 4 realises a 
chiral p x + ip y superconductor, although definitive evidence is still elusive [19] . An ultracold 
gas with tuneable p -wave interactions could be employed to explore the evolution from weak- 
to strong-coupling superfluids, including the topological quantum phase transition predicted 
in two-dimensional samples [20j ]2T] . In certain conditions, a p-wave superfluid could host 
Majorana modes that exhibit non-Abelian statistics, which is important for topological 
quantum computing mm- 

Here we report on time-resolved spectroscopic characterisation of an ultracold spin- 
polarised Fermi gas with near-resonant p-wave interactions. We find that both the momen- 
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turn distribution and the radio-frequency (rf) response follow asymptotic scaling consistent 
with p-wave contact relations [23H2S3, Furthermore, contact values measured with these two 
independent methods show good agreement. Then, relying upon the validity of the com¬ 
plete contact theory, we interpret the dynamical evolution of the contact as the population 
dynamics in the closed channel. 

Near-resonant p-wave interactions 

For p- wave scattering, the van der Waals interaction between atoms can be regarded as 
short-ranged [27]: i.e., acting only when r < r 0 , where r is the inter-nuclear separation and 
r 0 is comparable to the van der Waals length, about 3 nm for potassium. In ultracold Fermi 
gases, the inter-particle separation is much larger: 1/kp is typically 10 2 nm, where hkp is 
the Fermi momentum. As a result, the low-energy interaction can be described completely 
by the p- wave scattering phase-shift 5i m , fc 3 cot = —l/v m — k 2 / R m , where k is the 

relative wave vector, v m is the scattering volume, R m is the effective range, and m = {x, y, z} 
labels the projection of angular momentum onto the magnetic-field axis [271 [28] . 

Figure [I] depicts the control introduced by a Feshbach resonance [2]5], where a closed- 
channel molecular bound state i[) c is tuned near energetic resonance with an open-channel 
scattering state 'ipo of a pair of free atoms. Tuning is accomplished with a magnetic held 
due to a differential magnetic moment 8p between the closed and open channels. Far from 
resonance, v m and R m take on background values and R^, but near resonance, the 
scattering volume is resonantly enhanced: v m ~ m /(-£>o,m — B ), where B 0 rn and A m 

are the position and width of the resonance (see Methods and Fig. [I]). In contrast, R m > 0 
is only weakly field-dependent [231 ESI SO]- Since the resonance for m = z is split from the 
degenerate m = x and m = y resonances, both the strength and the anisotropy of p-wave 
interactions can be controlled with the magnetic held B = B z [16, 28}, [3T, 32] , 

For v m > 0, there is a Feshbach dimer state, which is a superposition of ^ c and i/j 0 , at 
energy E (im = —h 2 R m /Mv m , where M is the atomic mass. For v m < 0, the dimer state rises 
above threshold (E^ m > 0) and decays at a rate = 2 \[MR m Ef/^/h 2 (see Supplementary 
Information). At low energy, all p- wave resonances are narrow (fry m <$; Ep, where Ep is the 
Fermi energy) since 7 m /A'd,m decreases with energy. In some ways, they resemble narrow 
s-wave resonances [ 331 jMJ, which also have a quasi-bound dimer state above threshold. For 
the particular p-wave resonance we use, dipolar relaxation of -0 C to more deeply bound states 
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FIG. 1. Tuning p-wave interactions in 40 K. In a trapped gas of fermionic 40 K, the nature 
of interactions can be controlled by selecting internal states and tuning the magnetic field. Near 
the p -wave resonance studied here, spin-triplet collisions are resonant for atoms in the |2) state. 
This can be contrasted with the s-wave interactions near a Feshbach resonance at 202.2 G, where 
spin-singlet collisions in a |1) + \2) mixture are enhanced. Unlike s- waves, p -waves are normally 
suppressed at low energy due to a centrifugal energy barrier. The scattering volume v m diverges 
at distinct m = z and m = x,y Feshbach resonances, while R m is only weakly dependent on B 
and nearly isotropic. A parameterization of v m and R m is given in Methods. 
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limits the lifetime of the Feshbach dimer to < 8 ms | 28l [35 |. 

Universal p-wave contact relations 

That p -wave interactions might be described with an £ = 1 analogue of the contact 
was conjectured in ref J23j and recently followed by a full theory in refs EH and E3 The 
most significant structural difference from the s-wave contact theory is that each scattering 
channel has two contacts, U„, m and C R)rn , which are the thermodynamic “forces” conjugate 
to v and f?” 1 respectively. We summarize here what can be learned from the contacts, 
once they are measured or calculated. 

A. Thermodynamic identity. The change in free energy F for a uniform system is 

dF — —SdT + pdN - PdV-T j Y^C„, m dvP-T j Y^C Ft _ m dR- m \ (1) 


where S is entropy, T is temperature, P is pressure, V is volume, N is the total atom 
number, and /i is chemical potential. Note that C v _ rn has units of length, C Rtm has units of 
inverse length, and that both are extensive variables. In a harmonic trap, — P dV is replaced 
by 2(jJ~Ts C (U) du) osc , where (U) is the potential energy of the cloud and u> osc is the trapping 
frequency. 

B. Correlations. At short range, r -C kp 1 , the many-body wave function has a form that 
is controlled by two-body physics, but a normalisation that is controlled by the contacts [5] • 
For instance, the pair correlation function is 


g (2 \ r) -»■ 


'Cv kp C R)m /kp 
( k F r ) 4 ( k F r ) 2 


( 2 ) 


in the regime r 0 -C r kp 1 , where r = iq — r 2 is the relative coordinate, f = r/r, and Y\ m 
are the spherical harmonics for l = 1. 

C. Momentum distribution. The contacts also constrain the asymptotic form of the 
momentum distribution, n^. Indeed the s-wave contact is often defined to be the high-fc 
limit of k 4 nk [H Ej- For p- waves, in the asymptotic regime k R <C k l/r 0 , nj~ has two 
components [23ti26j : 

'v,m 2C R} m / q \ 

where k = k/k and V is the volume of the system. 


n k —> 


167T' 


V 




L(k) 
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D. Spectral response. The spectral weight of excitations that probe the high-energy or 
short-range sector of the many-body wave function are also controlled by the contacts j36l 
3D|. For rf transfer to a non-interacting probe state, the high-frequency tail of the spectral 
density is 



where ui is the detuning of the probe frequency from resonance, C v = J2 m C Vtm , and Cr = 



E. Fraction of the closed-channel molecules f c , m - Close to the Feshbach resonance, where 
v m » u„ g , fc,m is proportional to C v , m : 



(5) 


where i c ^ m = MS /i A m / h 2 p4j. In this aspect, C v _ m is similar to the s-wave contact [HE]- 
In contrast, Cn, m is an energy-weighted quantity that in the two-channel model also involves 
atom-dimer interactions (see Supplementary Information). 

Observation of the p- wave contacts 

The primary impediment to the exploration of p -wave many-body physics in trapped 
quantum gases has been atom loss that is faster or comparable to trap-wide equilibration [[551 
3M3]. Our experimental approach is to study the gas after a “quench” that quickly initiates 
enhanced p-wave interactions, accomplished with rf pulses. Before each pulse sequence, 
40 K atoms are confined in a crossed-beam optical dipole trap, spin-polarised in the lowest 
hyperhne-Zeeman state |1) and cooled to T « 250 nK, which is ~ 0.2 Af//cb, above the 
superfluid critical temperature [23) 31; 35j. A uniform magnetic held is stabilised at B = 
B 0m + dB m , in the vicinity of a p-wave Feshbach resonance for state |2). A resonant 40-ps 
7r-pulse transfers all atoms to 12), initiating tuneable p-wave interactions. After a variable 
hold time t, the gas is characterised either with rf spectroscopy or with time-of-flight (TOF) 
imaging, allowing contacts to be measured through relations (J4]) or ([3]) respectively, as shown 
in Fig. [2j Losses restrict t to be short compared to thermalisation times of low-energy or 
long-wavelength degrees of freedom. However, we find that spectra reach a quasi-steady- 
state, which likely reflects a local equilibrium. 

Radio-frequency spectroscopy probes the gas by transferring a fraction of atoms in |2) 


6 





FIG. 2. Observation of the p- wave contacts, a, The rf transfer rate T versus frequency uj at 
5B xy = +0.10(2) G and t = 160 ps shows a pronounced tail at u > 1 that can be fit to determine 
C v and Cr. Each point here corresponds to three repetitions of the time sequence shown as an 
inset. The transfer fraction is determined from an absorption image after Stern-Gerlach separation 
of states. b,c At the same 5B xy and t, the momentum distribution n(n) also has a visible k > 1 
(i.e. k > kp) tail, that is not present for a non-interacting gas (red line in c). The dashed line 
corresponds to the best-fit asymptote, C v ^ m kp/N = 0.028(6) and Cji_ m = 0.00. The distribution 
shown is the average of forty images, d, A comparison of contacts determined by T(cD) and by n{n). 
Circular markers are C V)Tn kp/N and square markers are Cfi_ m /kpN: blue markers have m = xy 
and red markers have m = z. Values determined at 5B xy = +0.10 are additionally outlined in 
black. The shaded region shows 1 s.d. uncertainty of the best-fit line, 0.96(7), and the dashed line 
has slope of 1. Error bars in a-d are statistical. 
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to the third-lowest energy state |3), which (like |1)) does not have resonantly enhanced 
interactions. The fractional transfer to 13), N 3 /(N 2 + N 3 ) = N 3 /N, is measured by state- 
selective absorption imaging, after a magnetic held jump that dissociates any Feshbach 
dimers. Figure^ shows an rf spectrum taken at 5B xy = +0.10(2) G. The transfer to |3) is 
given as a rescaled rate r(c5) = [Ep/h){Tikl 2 t Y: i)^ l {N 3 /N) 1 where cu = hu/Ep is the probe 
frequency rescaled by Ep, hVt/2 is the transition matrix element and Vtt r { is the pulse area. 
The latter is chosen to be small enough to probe the transition in the linear regime for 
u 1, where r(c5) oc /( c5). The high-frequency tail hts well to equation (Jr]) , and is used to 
determine C v and Cr (see Methods and SI for details). 

The momentum distribution is measured by resonant absorption imaging of the cloud 
after a 5.5 ms time-of-flight expansion. For these measurements, no rf pulse is applied 
at t ; instead, the held is rapidly jumped away from the p-wave resonance, preserving the 
interacting momentum distribution, which determines the ballistic flight after release from 
the trap. Figure^ shows the normalized distribution n(n) observed at 5B xy = +0.10(2) G 
versus k = \J {k% + /c 2 )//+, after azimuthal averaging in the image plane. Inherent to 
imaging is also a line-of-sight integration of n(k), so that the high-momentum scaling of 
n{n) is oc k~ 1 for the C V)m term and oc k+ 3 for the Cr^ term (Methods). Figure]^ shows 
that the leading order appears as an asymptotic plateau in k x n(n). A full fit with both 
terms is used to determine C v and Cr. 

Figure [iji compares the p-wave contacts determined from r(c5) and n(ft), across a range of 
magnetic held values. The dimensionless C v kp/N and CR/kpN are scaled by kp calculated 
using the peak density of a non-interacting gas, but results should be understood as an 
average over an inhomogeneous trapped ensemble pf6j. Since the contacts are only revealed in 
the asymptotic part of the distribution, analysis involves a low-energy cutoff, the systematic 
effect of which is studied in SI. Our analysis also assumes both C v and Cr are nonnegative, 
but the possibility of Cr < 0 is also discussed in the SI. 

The correlation between the two observables is 0.96(7), as determined by the slope of a 
best-ht line with no offset. This agreement, in addition to the observation of the predicted 
asymptotic scaling of equations (|3]) and Q, is strong evidence that the p- wave contact 
relations are valid. 

Field dependence of the p- wave contacts 

Figures |3^,,b show the p- wave contacts versus SB m near both the m = xy and the m = z 
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FIG. 3. The p-wave contacts near two Feshbach resonances, a, C v kp/N and b, Cft/kpN, 
versus magnetic field 5B m (lower axes) or dimer energy E ( \_ m J Ep (upper axes). Data are shown 
from both the m = xy (blue) and m = z (red) resonances, from rf spectra (filled) and momentum 
distributions (open). Momentum spectroscopy is shown for a smaller range of 5B m , due to limited 
signal-to-noise. Most data is taken at t = 160 /us, however long-time asymptotes from fits to 
dynamical data (see Fig. |4j) are also shown as smaller filled points, c, Numerical integration of 
measured C v gives the shift of free energy A F = F — F^g due to near-resonant interactions (see 
Methods for the details of the numerical integration). Data is referenced to Fy )g = F(B mSLX ) for 
5B m > 0; and to -^bg = F(B m i n ) for 5B m < 0. Illustrations depict the dimer energy, compared to 
the range of collision energies available in the Fermi sea. Error bars are statistical; see text for a 
discussion of systematic uncertainty. 
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resonances. The data includes contacts determined at t — 160/rs from F, at t — 160/rs from 
n, and asymptotic values from T versus t. The variable-t data (discussed in more detail 
below) also identifies a loss-dominated regime for 0.00(2) G < 5B m < 0.04(2) G, outside of 
which contacts reach a steady-state value despite atom loss of up to 20%. 

We observe a pronounced asymmetry about each Feshbach resonance: significant contacts 
are only observed for SB m > 0. C v is largest close to resonance, decreases with SB m , and 
vanishes beyond 5B m ~ 0.3 G, where E^m/E-p ~ 2. Cr instead peaks at SB m ~ 0.3 G before 
abruptly falling to zero for larger fields. 

Some of these salient features can be explained by a simple model, in which A T d = 
Y^ rn Nd,m non-interacting closed-channel dimers are in equilibrium with Nf free fermions. 
Each dimer has C v , m = 2 R m and Cr jm = —2 but free fermions make no contri¬ 
bution to the contacts. Since the m = xy and m = z resonances are well separated, 
C v kp/N ~ 2kpR m (N d /N) and C'R/{kpN ) ~ 2kpR m (Ed/Ep)(Nd/N). The assumption of 
equilibrium gives N d = (iV/2)(l — (E djin /2Ep) 3 ) in a harmonic trap at zero temperature [45 ] . 
This model would predict that both C v and Cr are the same near the xy and the z reso¬ 
nances, that C v — y 0 and Cr —>■ 0 as E dm — > 2Ep, and that a fully dimerised gas would have 
C v kp/N ~ 0.04, since kpR m ~ 0.04 in typical conditions. The additional factor of (E^m/Ep) 
in Cr gives Cr = 0 at resonance and at peak C'R/{kpN) ~ 0.06 at E d ,m/Ep ~ 1.6. 

Although this model does explain the peak value of C v and the range of 5B m at which 
significant contacts are seen, it does not explain the value or location of the maximum in Cr. 
A more realistic model would include finite temperature, and interactions between dimers, 
between atoms, and/or between atoms and dimers. For instance, resonant enhancement of 
atom-dimer interactions have been seen in a three-body calculation [SUED]. 

Independent of any particular microscopic model, but assuming adiabaticity, we can 
understand the thermodynamic implications of the observed contacts using equation (JTj). 
The change in free energy F versus 5B m is given by the integral of C v over v m x , assuming 
all other variables are constant. The contribution of Cr is not significant (Methods). The 
inferred A F is shown in Fig. [3]:-. The values shown have several possible systematic errors. 
First, some of the other variables that determine F are varied by 5B m : N decreases due to 
loss, and T increases by ~ 0.05Ep/kp near resonance. A second and more significant error 
may lie in the calibration of number and rf power, which combine to give a 30% systematic 
uncertainty in A F = F Fkg. Finally, equilibration is likely to be only local, and not trap- 
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wide. Despite these uncertainties, the integrated data is sufficient to demonstrate several 
qualitative regimes: 

(i) Below resonance (E d < 0), the gas is weakly repulsive, with 0 < A F <C E F . Here, 
resonant scattering is inaccessible to free particles, and the gas remains on the “upper 
branch” 13308]. Few or no dimers are formed, because energy-conserving two-body collisions 
cannot produce a dimer with a finite binding energy. Instead, the gas has weakly repulsive 
p-wave interactions. 

(ii) At resonance, we do not extract a value for F, because a steady-state in C v kp/N 
is not achieved, as discussed in the next section. However, the discontinuity in F between 
regime (i) and regime (iii) implies that the systems shifts from upper to lower branch in this 
region. 

(iii) Above resonance, in the range 0.25 < E d /Ep < 2, we infer a reduction in F per 
particle approaching half the Fermi energy near resonance. The significant reduction of 
free energy is partially explained by the formation of dimers, whose binding energy could 
contribute up to A F = —3Ep/A in a harmonic trap [35]. Additional contributions to 
A F include dimer-dimer or atom-dimer interactions. Accompanying the large A F in this 
regime are the largest observed contacts, and therefore the strongest p-wave correlations, as 
described be equation (|2]) . 

(iv) Farther above resonance, the scattering resonance at E d exceeds the maximum colli¬ 
sion energy 2 Ep in a zero-temperature Fermi sea, leaving primarily non-resonant interactions 
between atoms. In this regime, p-wave interactions are weakly attractive: 0 < —A F <C Ep. 

Dynamics of C v after the quench 

Evolution of the many-body wave function is required to adjust to the p-wave interactions 
initiated by transferring atoms from |1) into |2). We are able to observe these dynamics in 
the range 0.02 G < SB m < 0.10 G, by varying the time t at which we measure the contacts. 
For an isolated pair of atoms, the lifetime of the quasi-bound state q" 1 would set the key 
time scale (between 0.2ms and 1ms here). However, we find that the contact develops at a 
rate twenty to thirty times larger than q m . 

From Fig. [i]we know that only C v contributes to T for 6B m < 0.1 G. Measurement of T 
at a single rf frequency (typically oj ~ 5) is used to determine C v kp/N. We also record the 
atom number N at each t. Figures [4^-f show the time-resolved measurements of C v kp/N 
and N at various 5B Z . We observe that C v kp/N rises to an apparent steady-state value, 
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SB = 0.04 G 0.05 G 0.10 G 




f( ms) f(ms) f( ms) 
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FIG. 4. Dynamics of the p- wave contact and atom number close to resonance. C v kp/N 
(a-c) and N(t)/N(0 ) (d-f) are shown for B = 198.83 G (a, d), 198.845 G (b, e), and 198.89 G (c, 
f). Error bars in a-f represent the standard deviation over six separate measurements. Solid lines 
represent best fits to equation ([b]) ; dashed lines consider instead a dimer-dimer loss term discussed 
in the Supplementary Information. N(t) is normalised to the extrapolated N( 0). We vary A m , Lfd, 
and C v /Nd to fit the solution of equation (Jb]) to our data and iteratively determine a self-consistent 
.Bo,m for either resonance (see Supplementary Information). We find Bo,z = 198.792(13) G and 
Bo,xy = 198.301(11) G, which are used to determine the 5B m values reported in this work, g, The 
dimer association rates A xy (blue circles) and A z (red circles) versus magnetic field. The open 
circles indicate a loss-dominated solution. Solid lines in g are fits to A xy and A z with assumed 
proportionality to y m , finding A xy l^ xy = 16.3(3.7) and A z l^ z = 9.9(2.5). Dotted lines represent 
the theoretical prediction of \Q^ xy (red) and 87 ^ (blue), respectively. Error bars in g correspond 
to a confidence intervals. 
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while N decays relatively slowly. The initial growth rate of C v kp/N increases with 5B m , 
until it can no longer be resolved, at 5B m > +0.10 G. 

Near resonance, the wave function of each Feshbach dimer is dominated by [0H1 5DJ, al¬ 
lowing us to interpret the contact dynamics through a multi-channel model and equation (|5]). 
Three closed channels of p -wave dimers are initially empty, but come to equilibrium with the 
initially populated open channel of atoms, but all contribute to the observed atom number 
N — Nf + 2 Nd,m- We calculate the time evolution of Nf and N dim with rate equations 
that omit coherence between the channels but include collisional loss: 


Nf — — (2A m Nf — 2q m Nd,m + Lid,m NfNd,rnj , 

m 

-^d,m f NfN& ,ra? 


( 6 ) 


where A m is the dimer association rate, and Lfd,m is the fermion-dimer loss coefficient. We 
constrain the model to have Nd, x — N d _ y = 0 near the z resonance, and N d , z = 0 near 
the xy resonance. We furthermore combine the degenerate modes at the xy resonance to 
N dxy = N djX + N dyi A xy = A x + A y , and j xy = ^ xyy . We fit the solution of equation ([6]) to 
the data, hireling good agreement (see Figs. [4^-f). The rapid growth of iV d m for 5B m > 0 
was observed directly in refs 02] and [51] 

With Lf d ,m = 0, equation ([6]) would lead to a dimer population that asymptotically tends 
towards N d , xy = (Aj • y /'y xy )Nf and N d , z = (A z /'y z )Nf. The associated equilibration rates are 


1 / Txy.z 'Joci/jZ T 2A X y,z- 


(7) 


With loss present, we can distinguish between regimes (ii) and (iii). If r m Lf d ,m •C 1, the 
system reaches a quasi-steady state from which it decays slowly. This scenario is realised in 
Fig. [4j3,e and Fig. [4]c,f, and corresponds to regime (iii) in Fig. [3jr. If, however, r m L{ d/in 1, 
the loss rate of the system is faster than the association rate, and equilibration between 
N d and N { is inhibited. As a consequence, N d reaches a significantly reduced steady-state 
value on a time scale that is dominated by the loss rate. In turn, atom loss occurs on the 
time-scale of the dimer association. This case, corresponding to regime (ii) in Fig. |3ji, occurs 
closer to resonance (see Fig. [4^, d and data at 5B Z = + 0.04 G and SB Z = + 0.03 G in Fig. [ 4 ^) , 
because 1 /r xy , z — » 0 for 5B m —$■ 0 . The SI provides further details of how these two regimes 
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are distinguished, and demonstrates robustness against details of the loss model. 

Figure [3^ includes the asymptotic values of C v (smaller points). They are within un¬ 
certainty of C v at t — 160/us for SB m > 0.1 G, but provide an upwards correction for 
5B m ~ 0.05 G. With the inclusion of this correction, C v decreases monotonically with 5B m 
throughout regime (iii). 

Using only data from regime (iii), we find that A m are proportional to y m , with best-fit 
ratios A xy /2 r y xy = 8 .1(1.8) and A z /= 9.9(2.5) (see lines in Fig. |4jg) . These ratios are 
consistent with a perturbative treatment of resonant closed-channel molecular formation in 
a T = 0 Fermi cloud that predicts A m —>■ 87 ™ as 1/v —> 0 - (SI). The consistency between 
the dynamical response of C v and a model of dimer population supports the validity of 
equation ([5]). 

Conclusions 

While studied here for a Fermi gas in a metallic state, the p -wave contacts are expected 
to be “universal” in that they hold for any type of particle, whether boson or fermion, 
in any dimensionality, and in any state (superfluid or normal), so long as interactions are 
short-range and p- wave. Further tests of universality should include comparison to direct 
measurements of energy, structure factor, and dimer number, as has been done for the s-wave 
contact Hi 1521 fool. 

We have also identified a regime where the atom-dimer equilibration is much faster than 
loss. After equilibration, strong p-wave correlations persist for at least half a millisecond, 
eventually limited by dipolar and three-body loss rates. Searching for the onset of pair 
condensation in this dynamical window at a low temperature would be of great interest, 
especially in two-dimensional systems, where some loss mechanisms are less significant [50], 
and the possibility of chiral order exists [ 20 , 21 ]• 
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Methods 

Sample preparation. Spin-polarised 40 K atoms are cooled sympathetically with bosonic 
8 'Rb atoms. The |1), 12), and |3) states of 40 K refer to the high-field states adiabatically 
connected to the low-field rrif = —9/2, —7/2, and —5/2 states of the / = 9/2 hyperfine 
manifold of the electronic ground state, where / and rrif denote the total angular momen¬ 
tum and the corresponding magnetic quantum number respectively. At the end of cooling, 
residual 8 'Rb atoms are removed with a resonant light pulse leaving only 40 K atoms in state 
11), held in a cross-beam optical dipole trap. For spectroscopic data, N = 3.6(4) x 10 4 , 
while dynamics data was taken with N = 2.5(5) x 10 4 . The 7r pulse initialising dynamics is 
~ 98% efficient, so that we assume N 2 — N\. The mean oscillation frequency in the trap is 
uj osc /2tt = 440 Hz. 

The ideal-gas Fermi energy in a harmonic trap is Ep = hu osc (6Ny/ 3 , here ~ h x 30 kHz, 
and is also the local Fermi energy at the centre of the trap. Thermometry is based on 
fitting an ideal-gas Fermi distribution to an absorption image taken at high magnetic field 
after release from the trap. To measure the temperature after a 160/is hold time in |2), 
we jump the field rapidly to 209 G, turn off the trap, and transfer all atoms to |1) during 
time-of-flight. At 8B xy = 0.20 G, the apparent temperature rise is ~ 50 nK. Accompanied 
with a 25% number loss, this is an increase in reduced temperature k-pT/Ep of ~0.05. We 
note that with such short hold times, the gas is likely out of thermal equilibrium. 

State-resolved absorption images are taken at high magnetic field, after Stern-Gerlach 
separation of states 11), |2), and |3). Subsequent state transfer in a residual gradient field 
is used to map any population of interest to 11), since neither of the other states have an 
accessible cycling transition. For spectroscopic measurements, we optimise and calibrate 
this sequence to determine N 2 and N 3 . 

The magnetic field values are calibrated with the frequency of the 11)—12) transition, 
measured spectroscopically with a 120-//S pulse to an accuracy of ~ 1 kHz. The short-term 
accuracy of field values is ~ 2mG. Over one spectroscopic or dynamic data series, the 
magnetic field drifts by less than 10 mG. 

Analysis of rf spectra. The high-frequency tail in the spectral density I(u) of the rf 
transition to a non-interacting probe state (fl(u)du = N ) is given by equation (JI]). We 
measure the transfer rate T(u;) = N 3 (u)/t r f, with i rf the length of the spectroscopy pulse. 
The transfer rate is proportional to the spectral density and obeys the sum rule f r(u) du = 
Q 2 7 tN/2. We therefore identify r(cu) = Q 2 ttI(oj) /2. The scaled frequency u and normalised 
transfer rate F as defined in the main text are normalised to f F(O) du = 1/2. 

With these definitions, the high-frequency tail of the normalised transfer rate is 


1 C v kp ——1/2 , ^ Cr ~_3/2 

2 3 / 2 7t N 2 3 / 2 7t kpN 


( 8 ) 


Close to the p -wave resonances where C%/fcf <C C v , we measure C v using N 3 (u)/N at a 
single frequency value u: 


Cykp 

_ 

N 



(9) 


This single-frequency measurement technique is used in the data presented in Fig. [4} 

We calibrate the strength of rf fields by driving Rabi oscillations between |1) and |2) at 
209 G. At the peak power used in this experiment, we find a Rabi frequency of approximately 
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70 kHz with a typical uncertainly of 15%. 

The spectroscopic pulse has a Blackman envelope to minimise the frequency sidebands 
in the spectra. For a given pulse length t T f, the Blackman pulse area is a factor of ^0.4266 
different from that of a square pulse with the same length and peak power. The effective 
Rabi frequency that is used to normalise our rf spectra is Q/2n = 30(5) kHz. Within the 
measurement window, we find JT(uj)duj ~ 0.41(0.44) for 5B xy = +0.05 G (-0.10 G), which 
agrees with the sum rule to within systematic error. Note that when we take normalisable 
spectra the Rabi frequency must be reduced for \uj\ < 1, to keep the coupled fraction small. 


Analysis of momentum distributions. The momentum distribution in Fig. [2]d,c is obtained 
from a time-of-flight absorption image. The imaging beam propagates along the z direction, 
parallel to the Feshbach held. The optical density of the cloud is proportional to the column- 
integrated density p(x, y, t TOF ) = f dzp(r, ^tof) • For long time of flight, the initial cloud 
size is rendered unimportant and p(r, £tof) oc n(k). We normalize p by its pixel sum to 
obtain h(k) with f n(k)d 3 k = 1. Signal-to-noise is further improved by azimuthal averaging, 
yielding a distribution n(n) , where k = + k y /kp. Thirty to forty images are taken, 

and averaged. The contacts are determined with a fit to 


ti(k) -+ 


3 C v kp 
Hr N 


+ 


3 Cr. z + 3 C R>xy 
87t Nk F 


hi 


-3 


( 10 ) 


While the coefficient of the leading order is proportional to C v>m , just as in rf spectra, the 
sub-leading order is not simply proportional to We interpret the data assuming 

that C R)Z = 0 near the xy resonance, and that C R . xy = 0 near the z resonance. 


p-wave scattering parameters. For the p-wave Feshbach resonance in state |2) near 198.5 G 
we parameterise v rn and R rn as a function of magnetic held magnitude B using 


= * (f * or 

r SB m f5B m y c (SB m y 
Vm ~ v h r% A m + V A m ) \& m ) 




1 + 


6B„ 


A 


+ o 


R,m 


5B n 


A 


R,m 


2 


(ii) 


( 12 ) 


where SB m = B — R 0 ,m- From refEHj we use the values = (101.6 a 0 ) 3 , = (96.74 a 0 ) 3 , 


A, = 


and 


21.95G, A xy = 24.99 G, R h / = 47.19 a 0 , = 46.22 a 0 , A RtZ = -18.71 G, 

A RtXy = —22.46 G. From our own measurements, we use B 0)Z = 198.792(13)G and B 0tXy = 
198.301(11)G. Equations (11) and (12) are matched to ref [28], neglecting the next order 


in SB m /A. However, since \5B m \ < 0.4 G in this work, the next-order term has a relative 
magnitude 10 -6 for v and 10~ 3 for R. The binding energies measured in ref [35] matches the 
calculated Avb g /R hg to 3%, but further measurements are needed to constrain, for instance, 
A r . Our measurements of B 0z and B 0xy are consistent, within error, to those found in 
ref [25;, Both results may share a systematic offset due to the differential polarisability of the 
open and closed channel. 


Numerical integration to find AF. As a function of magnetic held, the reversible change 
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in free energy is 


dF 

dB 


E 


a f d(i/v m ) , a f d(i/R m ) 


+ 


< 9 ( 1/0 dB d(l/R m ) dB 


Assuming the xy resonance is well isolated from the 2 resonance, all terms in this equation 
can be evaluated with equation ([TJ) , measurements of and Cr^ and the parameterisa- 
tions pH) and (12). The dominant contribution comes from the scattering volume: 


dF ~ h 2 C v 
dB ~ 2MvmA 


m ‘-±m L 


1-2 


& Bm , /r\ ^rri \ 3 


+ 0(p) 

L-\ r 


(13) 


The change in energy between two magnetic held values, B\ and B 2 , is given by F(B 2 ) — 


" B 2 dF 
Si dB 1 

calculate the numerical integrals 


F(B 1 ) = f B 2 j^jdB. To evaluate the shift in energy due to near-resonant interactions, we 


cB 


F(B ) ~F(B ri 


F(B) ~F(B max ) - 


'B, 


min 

B n 


1B 


h 2 C v (B') 
2Mv h %A 
h 2 C v {B') 
2Mv hg A 


dB' for B < Bn 


dB' for B > Bn 


(14) 

(15) 


where we have suppressed the m subscripts for clarity, but are assuming that the z and xy 
resonances are isolated. To estimate the statistical uncertainty, we repeat the integration 
with data normally distributed around the mean values of C v at each B. The reported 
error is the standard deviation of the resultant F(B). The next order in SB m /A given in 
equation (13) can easily be included, but \5B m /A\ is < 0.01 when C v is above the noise 


floor of spectroscopy, so the resultant shift is not significant. 
The 5-dependence of F due to the change in R is 


dF 

dB 


-h 2 C 


R 


T.v 


2MR^n Ar^ 


(16) 


Across the range in which we observe p-wave contacts, this contribution to A F is small. 

|A|, even if 


One can understand this by comparing equations (13) to (16): 


since 


|A 


R\ 


Cr/Ic-fN is comparable to C v kp/N , the ratio of contributions to dF/dB near resonance is 
roughly /c|u bg /5 bg , which is ~ 1CT 2 for typical experimental parameters. 
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SUPPLEMENTARY MATERIAL 


Connection between closed-channel molecules and p-wave contacts 

Let the creation operators for molecules and fermions at position R be 0t(R) and ijr (R), 
respectively, and the internal wave function of the molecules be given by g( r) = g(r)Y\ m (f), 
with r the relative coordinate. Then the interaction Hamiltonian that converts scattering 
fermions to molecules can be written as 

Hint ~ 9 f drdR^(R)g(r)^R+ r -)^(R- r -) + E.c. (SI) 

where g is the (un-renormalised) coupling constant. Transforming to Fourier space and 
noticing that g(r) is an odd function, which requires its Fourier transform to be linearly 
proportional to kY lm (k), the effective coupling Hamiltonian can be written as [3D [57, 58j 

Hint = Hy 'y ^ ^^lm(^)^p a p/2+k a p/2-k H- C - (^2) 

V k,p 

where (ajj creates a molecule (fermion) with momentum p, and g is the coupling constant. 
We shall relate the latter to the scattering volume v m and also the effective range R m , for the 
scattering channel characterised by the magnetic quantum number m. The free Hamiltonian, 
on the other hand, can be written as 

^o = E< k )^k + E^(p)^ ( S3 ) 

k p 

where E( p) = p 2 /(4 M) + u m (B), where is m (B) = 5p m {.B — F>o,m) + const, is the detuning 
of the closed-channel molecule. S/i m is the magnetic moment difference between the open- 
channel scattering state and the closed-channel p- wave dimer. 

The connection of g and v m to v m and R m can be established by calculating the fermionic 
scattering T-matrix. The two-body Schrodinger equation can be written as 

(E - H 0 )^ = H m pV, (S4) 

and can be solved conveniently as T = T 0 + (E — R 0 ) -1 -£hnt'I / = T 0 + (E — Ffo)~ 1 7 1 T 0 , where 
To is the asymptotic state in the absence of interaction and we have defined the scattering 
T-matrix, T m . Developing, as usual, the perturbative series for T m , one obtains, in the 
centre-of-mass frame 


1 

T m (E) 


E - Um(B) 


2 g 2 ^ k 2 \Y lm (k)\ 2 

V ^ E + iO+ -tf/M 1 


(S5) 


where we note that the factor of 2 arises from the exchange symmetry of the identical 
fermions. To evaluate the integral, we need to impose an ultra-violet cutoff A in momentum. 
By further matching the relation with 1 /T m (E = h 2 q 2 /M ) -l/(v m q 3 ) - l/(R m q) - i, we 
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find 


1 _ 8n 2 h 2 

^ ” ~ Vm W + 

1 _ 87T 2 h 4 2A 

-Rm M 2 # 2 + 7T ‘ 


2A 3 

luT’ 


(56) 

(57) 


These two equations establish the relation between the two parameters in the Hamiltonian 
{g, u m } and the physical parameters {v m ,R m }. 

To obtain the number of closed channel molecules, it is simplest to use the Hellmann- 
Feynman theorem, by taking the derivative with respect to u m . That is 


Nd ,m ~ 


dE 

dUm 


dE dvj 
dv~^ dv 

uu m L/ty m 


(S8) 


Using the thermodynamic relation for C v , derived from equation (|Tj) in the main text, and 
the fact that close to resonance, u" 1 = — (B — B 0jTn ) / (A m v^f), we find that the number of 
dimers A d . m can be written as pREj 


Ab = 

1 v d,m 


h 2 


2M6^A m v^ 


a, 


(S9) 


In fact, close to resonance, since R m ~ 1/A ~ r$ with ro being the size of the closed channel 
molecules, the coefficient is just of order 1 /R m - Namely, 


C V m ^ 2 R-m A7], m for ^m ^ SlOO 


(S10) 


such that C v , m is directly proportional to the number of p -wave dimers in the closed channel, 
and the proportionality constant is of the order of the effective range R m >0. As a result, 
monitoring the value of C VtUl through rf spectroscopy is equivalent to monitoring the number 
of dimers in the closed channel. 


Likewise, combining the thermodynamic relations with equations (S6) and (ST), we find 


h 2 c R , m 

2 M 


dE 


dR- 4 


dE 




dvm 
M 2 g 


S dRdn 1 
2^2 


+ 


dE 


dg 


-2 


dg 


-2 


dR-} 


Rm 


g n 2 fy 4 \ U mEd,m+ 


(Sll) 


where H int is given in equation (S2). While the first term is related to the number of closed 
channel molecules, the later term involves the expectation value 


kp a p/2+k a p/2-k 


(S12) 


which describes the atom pair-molecule coherence. 


The model combining equations (S2) and (S3) has been used to study p -wave superfluidity 
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FIG. SI. Lowest order self-energy diagram for the dimer-dimer propagator. The vertex 
is given by gkY\ m (k)/y/V. The double line is for dimer-dimer propagator and single lines are for 
fermion propagators. 


in Fermi gases |H E3 E 8 ]. At zero temperature, at the mean field level, if condensation of 
molecules occurs only in the mth partial wave, one can show Cn. m = 2 /j,C Vjm M/h 2 , where /j 
is the chemical potential conjugate to the total number of fermions. The difference between 
2 pC V}m M/h 2 and —2R 2 n N dtm /v m is the contribution to C^ m from (H int ) in this case. 

Decay of quasi-bound p-wave dimers 

Our next task is to investigate how a p -wave dimer above resonance (E d:Tn > 0) decays 
into scattering fermions. It is simplest to compute the self-energy of the dimer propagator 
T, bjm (iQ n , p). The lowest order diagram is given in Fig. 

Explicitly, when analytically continued to real frequency 0, 


SI 


£ fe)m (0-H0 + ,p) 


2 .9 2 J 2\ v /(£P/2+k) + /(£p/2—k) - 1 

V ^ 1 lmUI nn + i o+-£ p/2+k -£ p/2 _ k , 


(S13) 


where £ k = h 2 k 2 /2A4 — p, and f(x) is the Fermi distribution function, describing the effects 
of the Pauli exclusion principle. In the two-body case (vacuum scattering), /i = 0 and the 
imaginary part of the self-energy is given by 


h7 m (0, p) = 2ImE ft m (n + i0 + , p) 

Mg 2 kMQ _p 2 \ 3/2 
47 i 2 h? y h 4 ) 


(S14) 


For p = 0 and when hVt = E dm > 0, the imaginary part is non-zero and is given by 

Mg 2 (ME d , m \ 3/2 


— ^'7m(-^'d,m5 0) 


47T 2 ^ 2 


- 


K 2 


= 2 


h 2 R m fME dt7n \ 3/2 


M 




K 2 


(S15) 


In the main text this relation is given as 7 m = 2 \fMR m E 3 J 2 n /K 2 . 

Analysis of the spectroscopic data 

Selection of the high-frequency tail of the rf spectrum. Since the contact only describes 
the functional form of rf spectra in the frequency range cu 3> 1 , we analyse only data above 
some cut-off frequency co > A, typically A = 3.5. 

There are two magnetic field regimes in which the fit results are insensitive to A. For 
5B m < 0.1 G, the spectra are dominated by C v ; whereas for SB rn > 0.2 G, the spectra are 
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FIG. S2. Effect of cutoff in analysis of rf spectra. The p- wave contacts a, C v and b, Cr for 
various values of A plotted as a function of magnetic field. The data plotted here is analysed with 
A = 3.0,3.5,4.0. The error bars are determined in the fit for A = 3.5. 


dominated by Cr. However for values of SB rn in between these two limits, the relative weight 
of the oH 1 / 2 and uA 3//2 components of the spectra is sensitive to A. 

We report data with A = 3.5 in Fig. [3] of the main text where the error bars show 
statistical errors from the fit. We repeat our analysis for A = 3, 3.5, 4 as shown in Fig. [S2j 
There is some scatter in the fitted values of C v and Cr, although it is roughly the same 
magnitude as the error bars from the fit. The reduced scatter in the data at the m — z 
resonance is due to increased sampling of the spectra. The general behaviour of C v and Cr 
are, however, not affected by the values of the chosen cutoff. With improved signal at large 
detuning we could minimise this region that is sensitive to the choice of A. 

Determination of the asymptotic power law. At each field value we fit equation Q from 
the main text to the high-frequency tail of our spectroscopic measurements of F. The high- 
frequency cut-off is chosen as discussed above and the prefactor to each power law is left 
as a free parameter. If we choose only one power law, setting the other prefactor to zero, 
the fits in the regime 0.1 < SB < 0.2 are poorly constrained. Further evidence that both 
power laws are required to fully describe our data is motivated by Fig. [3] in the main text. 
There are two regimes in SB where either only C v or Cr contribute to the spectrum. For 
SB < 0.1 G is dominant and therefore only the cu -1 / 2 power law scaling is present in the 
spectrum. This is shown in Fig. S3i at SB Z = 0.05 G. Fitting the high-frequency tail of 
this spectrum with a function, /( uj) = Auj we find an exponent r/ = —0.43 ± 0.06. This 
power law is made more apparent when plotting the high frequency tail of the spectrum in 
a log-log plot. The power law appears as a constant slope shown as the shaded region in 
Fig. [S3} > and compared to the theoretical slope of —0.5 shown as the dashed line. 

In contrast, for SB > 0.2 G Cr is the dominant contribution to the spectrum. Ex¬ 
perimentally, this appears as a change in the fitted value of r/. At SB xy = 0.3 G we find 

Qualitatively, the difference in power 


rj = —1.35 ± 0.21 from the data shown in Fig. 
law can be seen when comparing Fig. [S3|i and c. The weight at high frequency in Fig. [S3 


goes to zero much quicker as uj is increased than in Fig. |S3p .. When plotted on a log-log 
plot (see Fig. S31) the high frequency behaviour is clearly different. The fitted power law 
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FIG. S3. Determination of the power law in two limits of (iB. a, At 5B Z = 0.05 G the 
spectrum is dominated by C v . b, A log-log plot of the high frequency tail. The shaded region 
shows the fitted power law —0.43 ± 0.06 and the dashed line shows the prediction of —0.5. c, At 
5B xy = 0.3 G the spectrum is dominated by Cr. d, When plotted on a log-log plot the fitted power 
law with slope —1.35 ±0.21 (shaded region) is consistent with the prediction of —1.5 (dashed line). 

(shaded region) appears as a slope consistent with the theoretical value of —1.5 (dashed 
line), s-wave interactions would show a similar spectrum, but could only appear due to an 
unintentional admixture of incoherent atoms in state |1). We estimate this background to 
be less than 2% of the observed signal. 

Analysis of the momentum distribution data 

To improve the signal-to-noise in this measurement we perform several forms of averag¬ 
ing. Firstly, we average approximately 40 line-of-sight images. To perform this average we 
determine the centre of each atom cloud through a two-dimensional gaussian fit. We then 
overlay these centres to construct an average image. Secondly, we perform a radial aver¬ 
age. This serves the dual purpose of improving signal and averaging out noise. Using the 
centre determined in the previous averaging, we calculate the radial distance to each pixel, 
n = We then calculate the average of all pixels at the same radius. Finally, 
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FIG. S4. Effect of low-momentum cutoff in analysis. The p -wave contacts C v kp/N (circles) 
and Cu/kpN (squares) extracted from the momentum distribution at 5B xy = ±0.15 G as a function 
of momentum cut-off. A plateau in the fitted values appears around k c ~ 1.2. 
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we bin pixels to construct the distribution shown in Fig. [2jo, which has sufficient signal to 
perform a two parameter fit at large momentum. 

At each magnetic held the momentum distribution measured is highly dependent on the 
atom number which varies with the detuning from resonance (as the loss rate changes). As 
such, we independently choose a momentum cut-off k c for each data set. To determine the 
momentum cut-off we perform a power law fit using equation (J3| from the main text and 
vary the value of k c . The result is shown in Fig. S4 for data at 8B xy = 0.15 G. As can be 
seen, a plateau occurs over a range of cut-off values near k c ~ 1.2. We typically choose the 
mean value of this plateau as the cut-off for a given data set. 

Systematic effect of allowing negative Cr 

The values returned in the fits of our data with equations ([3]) and Q from the main 
text are sensitive to noise, atom loss at high energy, and small offsets due to imperfect 
background subtraction. These systematics combine to provide an uncertainty in the values 
of Cr, especially in the regions dominated by a large C v . Figure [So] shows values extracted 
for the contacts with Cr unconstrained. We find small negative values for Cr near resonance 
that are not consistent with zero. 

Our constrained fits in the main text are motivated in part by the fact that a small offset 
added to the spectrum can result in values of Cr consistent with zero. This offset is typically 
~ 50 atoms, consistent with the scatter in our number counting when no atoms are present. 
Further motivation for performing constrained fits comes from repeating this analysis on 
high-frequency tails simulated with Gaussian noise matching experimental conditions. For 
the maximum values of C v observed in this experiment we cannot statistically distinguish 
between a negative value of Cr and a small positive offset. 

We compare the results of these unconstrained fits for data from both the momentum 
distribution and spectroscopy in Fig. S6 A linear fit constrained to pass through the origin 


finds a slope of 0.86(9) (grey area), where uncertainty is statistical; compared to 0.96(7) in 
the main text. 
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FIG. S5. The p-wave contacts allowing negative Cr. a, C v and b, Cr for various values of 
5B with C v > 0 and Cr unconstrained. Negative values of Cr are observed near resonance. 
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FIG. S6. Comparison of the contacts, allowing negative Cr. A comparison of C v kp/N or 
Cn/kpN from both the momentum distribution and spectroscopy. Here, Cr < 0 is allowed, unlike 
in the main text. 


Rate equations near a p-wave resonance 

In the two-channel model described above, the dynamics of C v can be understood by 
considering how the dimer population changes as a function of time. Let the momentum 
distribution of the dimers be given by n m (p) and assume that the conversion rate between 
dimers and scattering fermions occurs at the same rate as that given in vacuum, y m , irrespec¬ 
tive of the asymptotic energy of the scattering fermions. This is a good approximation when 
the p-wave dimers are only slightly above the threshold E& m > 0, or very close to the p -wave 
resonance. In this case, due to energy and momentum conservation, if the centre-of-mass 
momentum of the dimer is p, then the incoming scattering fermion must have momentum 
p/2 + k and p/2 — k, with |k| = ko and tPh^/M = —Ed, m . The magnitude of k is fixed 
by the dimer binding energy, while its direction can be arbitrary. In the vacuum case, the 
angular average gives rise to the decay rate of the p-wave dimer, y m . Thus we can write 
down the phenomenological rate equations for the populations of dimers and fermions, 


dn r 


IP) 


dt 


= In 


|k|=fc 0 


/p/2+k/p/2-k(lfW ra (p)) (1 /p/2+k)(l /p/2-k ) n m{p) i (S16) 


where // means the angular average over k. The first term describes the conversion of two 
fermions into a dimer, while the second term describes the opposite process in which a dimer 
disassociates into two fermions, corrected with Fermi and Bose statistics. 

At the beginning of the dynamics we assume that there is a Fermi sea of Fermi momentum 
kp. When E^ m > 2 Ep (k 0 > kp) with Ep = kp/2M , the conversion into the the dimers 
from the fermions cannot happen, since energy and momentum conservation cannot be 
simultaneously maintained. On the other hand, in the limit ko/kp <C 1, i.e., that the dimer 
is only slightly above the threshold, we can write 


dn m ( p) 
dt 


In 


f p/2 (l + n m ( P)) - (1 - /p/2)Vn(p) 


7m[/p /2 - (1 - 2/ p/2 )n m (p)]. 


(sir) 


In the experiments, we are interested in the total number of dimers, which is related to C v>m 
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Summing over p on both side of equation (S17), one finds 


dt 


—7m y ^[/'p/2 (1 2/ p/2 )n m (p)] 

p 


7m^d,m T 7m ^ ^ /p/2 T ^ ^ /p/2^m(p)- 
P 


P 


(S18) 


The first term describes the vacuum decay of a dimer, proportional to N^m- The second 
term describes the conversion of two fermions into a dimer, which when summed over p, 
gives 8^ m N{ at zero temperature. The last term arises due to the Pauli principle, which 
inhibits the decay of dimers into fermions when the final states are already occupied and 
therefore reduces the apparent dimer decay rate. This term is more complicated to handle, 
and we combine it with the second term to be described by an empirical association rate 

A m - 


We can repeat the same steps for the density of scattering fermions and arrive at a set 
of coupled differential equations for m and Nf. 


dNd m 

dt 

dN { 

dt 


7m^d,m T A m N{ 
y] (27 m iVd, m — 2A m N{) 


(519) 

(520) 


which conserve the total number of fermions N — iVf + 2^ m A' r d m , and where the sum 
includes those closed channels relevant for the respective resonance. The above arguments 
indicate in our low temperature experiment A m = 0 when E^ m > 2 E F} and A m ~ 87 ™ when 
Ed >m is only slightly above the threshold. 


It is worth emphasizing that the scaling of the association term with Nf is a result of the 
degenerate open-channel Fermi sea. In the Boltzmann regime, a similar calculation predicts 
the intuitive Nf scaling for a two-body process: the association part A m N { in equations (S19) 
and (S20) is replaced by 


N? 

A m N{ —* 7 m y 


4vr K 2 
Mk B T 


3/2 


exp 


h 2 R n 


v m Mk B T 


(S21) 


The physical reason for this scaling is that the dimer association is the reverse process 
of dimer decay, which in vacuum would occur at y m . The ratios are determined by the 
phase space density of the fermions that pair into dimers. In a thermal cloud, the ratio is 
proportional to nA 3 , where A is the de Broglie wavelength and n is the local density. For 
a degenerate Fermi cloud, the initial rate saturates to 87 at 1/v = 0 and T = 0. Thus the 
observed magnitude of A m reflects the degenerate nature of the cloud throughout the fast 
dynamics. 


Empirical atom-loss term. We consider two different inelastic collisions that lead to a 
loss of atoms from the trap: collisions between two dimers (L dd ) and between a fermion and 
a dimer (L fd ). Including these processes with respective loss rates L, equations (S19|S20) 
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FIG. S7. Analysis of dynamics data, a, Typical xf ed map (here: B = 198.87 G, corresponding 
to 5B Z = 0.08 G) showing two local minima corresponding to the association-dominated (i) and 
the loss-dominated solutions (ii). The dashed line shows the condition r~Lf d = 1. b, Best fit 
values for the association rate at both resonances, and for both loss models. We also include the 
discarded local solutions for the p z resonance, which clearly violate the predicted proportionality 
to 7 m oc (5B^ 2 (small symbols). Solid and dashed lines are as in Fig. Ug in the main text. 


take the form 


Na,m = -7mWd,m + A m N { - L™ N<N d , m - ^ ] N d , m N d , m , 

m' 

N { = -2 ^ [A m N { - lm N d , m + \L { ™ ] A f iV d>m ) . 


(S22) 


(S23) 


Near the p z resonance, equations ( S22|S23 ) provide two differential equations that we use to 
model C v (t) oc iV dj ~(t) and N(t). Near the p xy resonance [m = x,y), we have three coupled 
differential equations for Nf, N djX and N d i2r These can be reduced to a pair of equations for 


= 7 „, = L™, and L™ = L 


■ (y) 


(*) _ r(y ) 
dd — -^dd ' 

so that all 


Nf and N djXy = N djX + N d>y if one assumes A x = A y , -f x = ^ y , 

Our three-dimensional geometry ensures cloud radii that exceed (u m ) -3 and R mj 
these conditions are well fulfilled. We can therefore use the same set of two equations for 
the p z and the p xy resonances, where we interpret A xy = A x + A y in the latter case. 

Analysis of the dynamics data. For each sampled held value, we fit the solution of equa¬ 
tions ( S22||S23 ) simultaneously to our dynamical measurements of C v (t)kp/N(t ) and N(t). 
In this, we choose only one of the loss terms, setting the other loss rate to zero; varying more 
than one loss parameter leads to poorly constrained fits for our data. We furthermore as¬ 
sume a resonance position, which we use to calculate through equation (S15). This leaves 


us with four fit parameters: the pre-factor a = C v /N d , the phenomenological association 
rate A m , the respective loss rate L and the initial atom number Nq = N(t — 0) which we 
need to include due to an unknown calibration factor in our state-selective imaging scheme. 
For a fixed grid of pairs (A m , L ), we calculate the reduced chi-squared value yf ed (A m , L) for 
the optimal values of a and N d . 

A typical distribution is shown in Fig. [S7f i for a held value near the p z resonance. We 
typically find two local minima with similar values of x(? ed for each set of dynamical data. 
One corresponds to the loss-dominated regime described in the main text, while the other 
would have the system reach a quasi-steady state. Figure [S7| ) shows the association rates 
A for either a pure fermion-dimer loss term (while L dd = 0), or a pure dimer-dimer loss 
term (while Lf d = 0). Both loss models give compatible results. The result of the Fermi’s 
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golden rule calculation, A m oc 7 m , is employed to discard one solution for each time-resolved 
measurement (small symbols in Fig. |ST|b) . The remaining best-fit values for A m outside 
the loss-dominated regime are then extrapolated to A m = 0 to determine a new resonance 
position Bq"‘\ This new held value is assumed for a next iteration until the fit routine has 
converged. 
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